Лабораторная работа №11.

Символьные вычисления в Maple

Бесконечные и конечные суммы и произведения

Вычислим конечную сумму sum((1+i)/(1+i^4),i = 1 .. 10),

> sum( (1+i)/(1+i^4), i=1..10 );

51508056727594732913722/40626648938819200088497

затем определим значение бесконечной суммы sum(1/(k^2),k = 1 .. infinity)

> sum( 1/k^2, k=1..infinity );

1/6*Pi^2

Аналогично действуем для нахождения конечных и бесконечных произведений, например product((i^2+3*i-11)/(i+3),i = 0 .. 10)

> product( ((i^2+3*i-11)/(i+3)), i=0..10 );

-7781706512657/40435200

Вычисление выражений

Зададим следующее выражение (41*x^2+x+1)^2*(2*x-1).

> expr1 := (41*x^2+x+1)^2*(2*x-1);

expr1 := (41*x^2+x+1)^2*(2*x-1)

Бывает полезно уметь раскрыть скобки :

> expand(expr1);

3362*x^5-1517*x^4+84*x^3-79*x^2-1

затем снова сгруппировать переменные

> expr2:=factor(%);

expr2 := (41*x^2+x+1)^2*(2*x-1)

Вы можете также упрощать выражения самого разного типа. Например, упростить cos(x)^5+sin(x)^4+2*cos(x)^2-2*sin(x)^2-cos(2*x)

> simplify( cos(x)^5 + sin(x)^4 + 2*cos(x)^2 - 2*sin(x)^2 - cos(2*x) );

cos(x)^5+cos(x)^4

Можно упрощать выражение при помощи команды normal , она позволяет упрощать дроби, сокращая их.

Например, упростим следующую дробь (x^3-y^3)/(x^2+x-y-y^2).

> normal( (x^3-y^3)/(x^2+x-y-y^2) );

(x^2+y*x+y^2)/(x+1+y)

Обратите внимание, что вы можете вычислить значение выражения, содержащее переменные в некоторой точке:

> eval(expr1 , x=1 );

1849

Дифференцирование и интегрирование

Зададим функцию proc (x) options operator, arrow; x*sin(a*x)+b*x^2 ....

> f := x -> x*sin(a*x) + b*x^2;

f := proc (x) options operator, arrow; x*sin(a*x)+b...

Вычислим производную данной функции diff(x*sin(a*x)+b*x^2,x)и запишем в переменную f_prime.

> f_prime := diff( f(x), x);

f_prime := sin(a*x)+x*cos(a*x)*a+2*b*x

Вы можете также посчитать вторую производную: diff(x*sin(a*x)+b*x^2,x,x)

> f_second := diff( f(x), x, x );

f_second := 2*cos(a*x)*a-x*sin(a*x)*a^2+2*b

Для проверки полученного значения вы можете взять неопределенный интеграл от полученного выражения производной

> int(f_prime, x);

-cos(a*x)/a+(cos(a*x)+a*x*sin(a*x))/a+b*x^2

> simplify(%);

x*(sin(a*x)+b*x)

Для вычисления определенного интеграла достаточно в оператор int добавить область изменения переменой.

Вычислим int({sin(a*x)+x*cos(a*x)*a+2*b*x},x = 1 .. 2)

> int( f_prime, x=1..2 );

2*sin(2*a)+3*b-sin(a)

Вычисление пределов

Вы можете вычислять как конечные, так и бесконечные пределы для широкого класса функций. Вы можете брать как левый, так и правый пределы.

Вычислим следующий предел: Limit((2*x+3)/(7*x+5),x = infinity).

> expr := (2*x+3)/(7*x+5):

> limit( expr, x=infinity );

2/7

Теперь продемонстрируем вычисление левых и правых пределов: limit(tan(x+Pi/2),x = 0).

> limit( tan(x+Pi/2), x=0, left );

infinity

> limit( tan(x+Pi/2), x=0, right );

-infinity

Maple может распознать, когда предел не определен limit(sqrt(1-cos(x))/x*sqrt(2),x = 0).

> limit( sqrt(1-cos(x)) / x*sqrt(2), x=0 );

undefined

Степенные ряды

Maple с легкостью раскладывает аналитические функции в ряды Тейлора до заданного порядка

> expr := sin(4*x)*cos(x):

> approx1 := series( expr, x=0, 6 );

approx1 := series(4*x-38/3*x^3+421/30*x^5+O(x^6),x,...

Тот же результат можно достичь, используя функцию taylor

> taylor( expr, x=0, 8 );

series(4*x-38/3*x^3+421/30*x^5-10039/1260*x^7+O(x^8...

Задания:

1.      Вычислить: .

2.      Вычислить предел: .

3.      Продифференцировать функцию: .

  1. Найти аналитически точки пересечения графиков функций:  и . Результат проверить графически.
  2. Исследовать функцию на наличие асимптот: . Построить график для проверки.
  3. Разложить на множители выражение
  4. Вычислить площадь области, ограниченной кривыми   Image1174.gif (1080 bytes). Построить графики.
  5. Найти наибольшее и наименьшее значение для функции на отрезке [0,4]
  6. Построить касательную и нормаль к графику функции Image530.gif (1073 bytes)в точке (1, 0).
  7. Определить интервалы выпуклости, вогнутости и точки перегиба функций
  8. Разложить функциюImage672.gif (970 bytes) в окрестностиImage671.gif (954 bytes).